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Abstract. The aim of this paper is the analysis of the fractional Poisson process where the state 
probabilities p k k (t), 4^0, arc governed by time-fractional equations of order < Vfa ^ 1 depending 
on the number k of events occurred up to time t. We are able to obtain cxplicitcly the Laplace transform 
of p k k (t) and various representations of state probabilities. We show that the Poisson process with 
intermediate waiting times depending on differs from that constructed from the fractional state 
equations (in the case Uk — for all fc, they coincide with the time-fractional Poisson process). 
We also introduce a different form of fractional state-dependent Poisson process as a weighted sum of 
homogeneous Poisson processes. Finally we consider the fractional birth process governed by equations 
with state-dependent fractionality. 



1. Introduction 

We first consider a state-dependent time- fractional Poisson process N(t), t ^ 0, whose state probabilities 
p u k k (t) = Pr{iV(£) = k} are governed by the following equations 

" ^pl" (t) = -Xpl" (t) + Ap^- 1 (t), k > 0, t > 0, v k e (0, 1], A > 0, 
1 ' ■ ' ' lv hfn , fl k = 0, 



k>l, 



where p^ h (t) = 0, if k e Z~\{0}. These equations are obtained by replacing, in the governing equations 
of the homogeneous Poisson process, the ordinary derivative with the Dzherbashyan-Caputo fractional 
derivative that is [10] 

, , A v , s { rl 1 r t(t-s) m - l '- 1 f ( - m) (s)ds, m-l<v<m, 

(L2) *r m -\gz, 

We remark that in (1.1), the order of the fractional derivatives depend on the number of events occurred 
up to time t. By definition we have that 

(1-3) ^P h * (*) = \\t - ^ (s) 6s, 0<u k <l. 

Hence the dependence of (t) on the past is twofold. On one side, the fractional derivative depends on 
the whole time span [0, t] through the weight function. On the other side the number of events occurred 
up to the time t modifies the power of the weight function. This means that the memory effect can play 
an increasing or decreasing role, in the case of a monotonical structure of the sequence of fractional orders 

For example, if Vj, decreases with k, the memory function tends to be constant and to give the same 
weight to the whole time span [0, t] . We notice that state-depending fractionality was considered in different 
contexts by Fedotov et al. [5]. 

For i/fc = v, for all k, the system (1.1) coincides with the one governing the classical fractional Pois- 
son process considered by Beghin and Orsingher [3], where the fractional derivative is meant in the 
Dzherbashyan-Caputo sense as in this case. Of course, if v k = 1, for all k, we retrieve the governing 
equation for the homogeneous Poisson process. A quite different approach to the fractional Poisson pro- 
cess was developed in Laskin [6], where the fractional derivative appearing in the equations governing the 
state probabilities coincides with the Riemann-Liouville derivative. 

By solving equation (1.1), we obtain that 

(1.4) e- at p u , k {t)At = —f^ , s>0. 

Jo n-= (^+A)' 
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The inversion of (1.4) is by no means a simple matter and we have been able to obtain an explicit result 
for Pq° and p" 1 in terms of generalized Mittag-Leffler functions defined as (see for example Saxena et al. 
[12]) 

(1.5) W = |^TPn)' ,,/>,meR+«lt 

We give also the distribution p 1 ^ (t) of the Poisson process with fractionality Vk depending on the number 
of events k, in terms of subordinators and their inverses (see formula (2.26) below). 

A part of our paper is devoted to the construction of a point process Af(t), t ^ 0, with intertime Uk 
between the fcth and (k + l)th event distributed as 

(1.6) Pr{U k >t} = E Vk , 1 (-At Vk ). 
The Laplace transform of the univariate distributions of Af(t), t > 0, is 

(1.7) 1 e - st Pr{AA(t) = fc} di = A ft - F ^— -, 

Jo n J=0 ( s J + A ) 

which slightly differs from (1.4). From this point of view the state-dependent fractional Poisson process 
differs from the time-fractional Poisson process because the approach based on the construction by means 
of i.i.d. inter-event times Uk and the one based on fractional equations (1.1), do not lead to the same 
one-dimensional distribution. We show that the probabilities Pk(t) = Pr{7V(i) = k} are solutions to the 
fractional integral equations 

(1.8) Pk(t)- Pk (o) = -\r k P k(t) + xj"*- i p*_i(t), 

where I" h is the Riemann-Liouville fractional integral 

(1.9) (/"*/) (t) = \\t - sr-\f(s)ds, Vk > 0. 

1 (Vk) Jo 

A third definition of the state-dependent fractional Poisson process, say N(t), with distribution 

(At) 3 ' 1 

d- 10 ) Pr{JV(t) = j} = ^ +00 (At)j 3 j > 

Zjj=o r(^j+i) E Vj:1 {\t) 

is introduced and analyzed in Section 3. The distribution 

{\tf 1 



(1.11) p r {7V„(i)=j} = 



investigated in Beghin and Orsingher [3], has been proved to be a weighted sum of Poisson distributions 
in Balakrishnan and Kozubowski [1] and Beghin and Macci [2]. 

Finally, we analyze the state-dependent nonlinear pure birth process with one initial progenitor, where 
the state probabilities p^ fc (t) satisfy the fractional equations 



(1.12) 



^Pl h (t) = -\kPl k (t) + \k-ip7-Wt), k>i,t>o,u k E (o, i], 



KMO) 



[0 k^2, 

The Laplace transform of the solution to (1.12) reads 



(1.13) 



r+aa /k—l 

e- st pl*(t)dt= [Y\\ 
Jo Vi=i 



The case where Vk = v, for all k, has been dealt with in Orsingher and Polito [8]. A similar and more 
general state-dependent fractional birth-death process was recently tackled by Fedotov et al. [5], where 
possible applications to chemotaxis are sketched. 



STATE-DEPENDENT FRACTIONAL POINT PROCESSES 



3 



1.1. Notation. For the sake of clarity we briefly summarize the notation used for the different point 
processes analyzed in the following sections. 

First we indicate with N(t), t ^ 0, the counting process associated with the variable-order difference- 
differential equations (1.1). In particular the state probabilities p k k (t) = Pr{iV(t) = k}, k ^ 0, represent 
the probability of being in state k at a fixed time t ^ 0. The point process constructed and studied in 
Section 3 by means of independent but non i.d. inter-arrival times is instead indicated by a calligraphic 
J\f(t). Both processes, when = v for all k ^ 0, reduce to the time- fractional Poisson process N v (t) treated 
for example in Beghin and Orsingher [3] . In the same article it is also considered the alternative definition 
for a fractional Poisson process characterized by the distribution (1.11) and denoted here by N v {t), t ^ 0. 
We refer to its direct generalization in a state-dependent sense as N(t) for which the distribution becomes 
that in (1.10). Lastly, the linear fractional pure birth process with state-dependent order of fractionality 
presented in the last section is simply indicated as Nu n {t), t ^ 0. 

2. The state-dependent fractional Poisson process 

We first consider a state-dependent time- fractional Poisson process N(t), t ^ 0, whose state probabilities 
p k k (t) = Px{N(t) = k} are governed by equations (1.1). We have the following result 

Theorem 2.1. The Laplace transform of the solution to the state- dependent time-fractional equations 

( ^' (t) = ~Xp\ k (t) + Xpl'r, 1 (t), k>0,t>0,v k e (0, 1], 



(2.1) 



reads 



P* fc (0) = 



1 k = 0, 
k>l, 



(2.2) P k k {s)=\ e- st pl"{t)At = 



~Vu / \ I — st V 

Pk (*) = e p k 
Jo 



where the fractional derivative appearing in (2.1) is in the sense of Dzherbashyan-Caputo. 

Proof. We can solve equation (2.1) by means of an iterative procedure, as follows. The equation related 
to k = 



SaPoW = -Apb°(0. *>0, n,e(0,l], 



[Po (°) = !- 

has solution Po°(t) = E^ (u i(—Xt v °), with Laplace transform 

(2.4) pg»(,)= r + V st ^(t)d*= 

Jo 



A + s"o 
where 

(2.5) ^A^^Az^L 

k—0 

is the Mittag-LefHer function. 
For k = 1, the equation 

r9fi s ffei 1 (*) = -AK 1 (i)+Apo°W. t> 0,^6(0,1], 

( } K(o) = o, 

has solution with Laplace transform 

r+oo 

(2.7) P?(s)=\ e- st p?(t)dt 

Jo 



Xs" - 1 1 
A + s' y o A + s"i ' 

By iterating this procedure, we arrive at 



(2.8) P k k {s)=\ e- st pl»(t)dt = 



f + OO 

~vu / \ I —st v 

Pk 00 = e p k 
Jo 



X k s v °- 1 



n-.o^ + A)' 

□ 
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Remark 2.2. A direct approach based on the inversion of the Laplace transform of (2.8) is clumsy and 
cumbersome. We give the explicit evaluation ofp^it). In this case, from (2.6), we have that 

P?(s) = f +0 °e- s V 1 (i)dt 
Jo 

As" - 1 



(2.9) 



A 2 + A(s"o + s"i) + s"o+^i 
As" ' 1 1 

1 ~<~ A 2 + s"0 + "'l 



As" " 1 2 



(-A(s' y ° + s^ 1 ))™ 



CO , ,x m m / 

Zj f\2 i -S i/o+i'i W+l Z-J I 



n ( A2 + a" *" 1 ) n v - 

TTie inversion of (2.9) involves the generalized Mittag-Leffler function, defined as (see for example 
Saxena et al. [12]) 

{-\t") k ?{m + k) 



(2.10) 



^ fc!r(^fc + /3)r(m) ' 



where v,f3,me K + . 

Indeed, we recall the following relation. 

r+co 



(2.11) 



f 

Jo 



vm—P 



(A + S v ) r 



v {n 



In view of (2.9) and (2.11), we arrive at 

(2.12) p?(t) = 

for f/ie cose ^i = i^o = ^, formula (2.12) becomes 
oo m / \ 

^m-vm+l V" 1 I m 1 .^(m+l) / \2,2i/n 



go m / \ 

m = r=0 \ / 



(2.13) p?(t) = ^ 

m — C 

00 

= s 

m — C 

00 

m — C 

00 

00 

= 2 

m — C 

00 



-l) m (At") r 



J 2i/,i/(m + l) 



i )+ i(-A 2 ^) 



-l)° 1 (AO m+1 2 m V ( m + r \ (-l) r (A 2 t 2i y 



^ V r I T(2vr + v{m + 1) + 1) 
-l) m (At") m+1 2 m ^ I ~~ ^ J ( A 2 t 2l/ ) r - f e w w~ 2vr ~ v ( m + 1 )- 1 



dui 



-i/(m + l) — 1 



00 



e u> 

* ^ a T" — 

„,-"(m+l)-l 

e 



(AM 2 2^i J, 
~ ; ; 2^i J Ha " (\H 2 »w- 2 » + l) m+1 



-(m + 1) | 2.2k -2knt 

r | (A t to ) 



dto 



I* , — f — 1 



At" iw 



S Mr 



2w-"At" 



27ri JHa i. W " + At") 2 

At" 

E ViV (—\t ), 



dw 



where we have used in the last equality the fact that 



(2.14) 



da; 2-ki ax 



lew , v I e % 
aw = t 

JHa W" -X 27TI J Ha {W v 



-dw, 
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and we have applied the contour-integral representation of the reciprocal of the Gamma function 

(2-15) _L = -L f e"«-dtt, 

r(as) 2tu J Ha 

where Ha stands for the Hankel contour (see Beghin and Orsmgher [3]). 

We notice that equation (2.13) gives the result obtained for the time-fractional Poisson process in Beghin 
and Orsingher [3] as expected. Moreover by considering that 

r+oo \+v \ n u—\ 

(2.16) 

Jo 

we retrieve, for the case v = vo = v\ that 



(2.17) 



p?(t) = — E u . v {-Xt v ), 
v 



that is the result obtained for the time-fractional Poisson process (see formula (2.11) of Beghin and Ors- 
mgher [3]). 

By applying formula (34) of Saxena et al. [12] it is possible to give an explicit expression for p^ k (t) , for 
any k ^ 2, in terms of cumbersome sums of generalized Mittag-Leffler functions. 

Remark 2.3. A different way to give a representation of the state probability in the state- dependent 
Poisson process is given by the following integral approach; starting from (2.8), we have 



(2.18) 



'(«) = f 
Jo 



t p" k k (t)dt = A* 



nJ-o(^+A) 



I" 

Jo 



dvjQ 



nj 

3 = 1 J ° 



— \wa \ —WAS 3 . 

e 3 Ae 3 dw. 



■Pt{W(x) > t) = i-\ h v {s,x)ds, 
ox i t 



For the following developments, it is useful to recall that the inverse process of a u-stable subordinator 
H"(t), t^O, namely C v \i), t 5j 0, is such that 

(2.19) Pr{£"(i) < x} = Pt{U"(x) > t}, x,t>0. 

Hence the relation between the law l v (x, t) of the process C v (t) and the law h v (x, t) of the process H"(t) is 
given by (see for example D'Ovidio et al. [4]) 

(2.20) l Ax ,t)= Pl{C J edX} =l 

dx ox 

or otherwise 

(2.21) Pi{H"(x) e dw} =\ Pr{£ v (t)edz}. 

Jt Jo 

Hence the density of the inverse process C v (t) reads 

(2.22) Pr{C(t) edx} = f Pr{H v (x) e dw}. 

dx J t 

Therefore the Laplace transform ofl v (x,i) is given by 

(2.23) Ux,s)= \ e- s %(x,t)dt= f e" st -^-[f + Pt{H"(x) e dw} dt 

Jo Jo dx Ut 

d p 00 r w 
-°- Pr{H"(x) e dw} e^'di 
Ax Jo Jo 

TOO 

(l-e- sw )Pr{n u (x) edw} 

Jo 



dx 
1 d 
s dx 



where we used the fact that 
(2.24) 



,(x,s) = f 
Jo 



h v (x, t)dt 



We also notice that the explicit form of the law of the inverse of the stable subordinator is known in terms 
of Wright functions [4J. Going back to equation (2.18) and in view of (2.23), we can write 

k 

_-Xw 



(2.25) pl»(s) = 



J e- Xw °dwoj e- st L (w ,t)dtj ([] A J 



rco ' 

e~ w idu)j e~ sx h v . [x, Wj)dx 
Jo , 
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poo poo pc 

dwoe Xw ° ■ ■ ■ \ dw k e~ Xwk I 
Jo Jo Jo 



e st l„ (w ,t)dt 



fc poo 

FT A e~ ax h Vj (x,vij)dx 
3=1 Jo 



Hence by inverting the Laplace transform we obtain 

poo poo poo 

(2.26) Pl k (t) = ^ k \ dw e~ Xwa \ dwie~ Xwi ■ ■ ■ dw k e~ XvJk [l vo (w, t) * h viz ... ,„ fc (wi, ■ ■ ■ , w k , t)] , 

Jo Jo Jo 

where the symbol * stands for the convolution of the law of the inverse stable subordinator l VQ and the 
distribution of the sum of k independent stable subordinators h vi , — ,v k (wi, •• • ,Wk,t). In other words 
l VQ (w,t) * hux,-- , Bj L (wi, • • • ,Wk,t) is the distribution of the r.v. 

fc 

(2.27) C°(t)+^]H Vi (t). 

3 = 1 

Remark 2.4. Another interesting characterization of the state-probabilities of the above process is given 
by the following observation. First, we recall that 



(2.28) 



(2.29) 



t t v - 1 E v , v {-\t v )dt 



p+o 

Jo 

p + OO 

e- st E v , 1 {-Xt v )dt = 

Jo 



1 

~ A + s"' 



A + s" ' 



Hence, from (2.8), u>e /ma! f/iaf 

A«= s "o-i 



(2.30) 



v7 (s) 



f + °°e- st ^ 0il (-Ar)dtl n [f 

Jo J , =1 LJo 



On the other hand, from (2.25), we have 



(2.31) pl"(s) = 



e st l UQ {ii)Q,t)&t 



poo rcc 

e- Xwo dw 
Jo Jo 

£°e- s *L (A,t)dt) ffl A £ 



& roo roc 

j=i Jo Jo 



e sx h Vj (x, Wj)dx 



°h Vj {x, X)d3 



which clearly coincides with (2.30). 

By inverting the Laplace transform, we obtain the following result 

(2.32) Pl k it) = E voA (-Xt"°) * At^-\E^(-A^). 

3=1 

We notice that, the last equation can be written in terms of Prabhakar operator, that is an integral operator 
involving a Mittag-Leffler function as kernel [11]. From equation (2.32) we have an integral representation, 
in explicit form given by 

(2.33) Pi 1 W= f E V0 , 1 (-\(t-s) V0 )E vl , vl (-\a vl )8 , ' 1 - 1 d8 

Jo 

P2 2 W = f dBi^,i(-A(t - si)" ) da28?- 1 E vl , ul (-\8?)(ai-a a ) va - 1 E vt , Ua (-\(8i- Ba ) va ) 
Jo Jo 



vl k {t)= f 

Jo 



= dsxE^xi-Xit-sxf ) ds 2 



Jo 



x P dsks?- 1 1 E Vk _ lt „ k _ 1 (—Xs'^ h ~ 1 )(sk-i - sur-'E^i-Xisk-i - s k y»). 
Jo 

In order to find the mean value of the distribution p k k {t), we multiply all the terms of (2.1) for k and 
sum over all the states so that 



(2.34) 
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In the case Vk = v, for all k, we have 

At" 



(2.35) 



d" 



E(iV„(t)) = A, 



whose solution is given by H(N w (t)) = r ^J +1 ^ ( see formula (2.7) in Beghin and Orsingher [3]). We notice 
that it is possible to find an interesting summation formula by using the Laplace transform in equation 
(2.34). Indeed we have 



(2.36) 

and recalling that 
(2.37) 

we find that 
(2.38) 



^ ks Vk p v k k (s) = As 1 , 



Pk («) 



As 



This summation formula is not trivial and we can check that it works for example in the special case v = Vk 
for all k. 



(2.39) 



fcA fc s 2 



£i ( s " + s " + A ^ti (s" + A) 



2 



fcA fc 



As 2 



s" + A 

As 2 " 
s" + A 

As 2 " 
s" + A 



7 

a 

d™ £± {s v + A)* 
d w 



dw s" + A — w 

s" + A 
(s" + A - w) 2 



= A 



Remark 2.5. VFe notice that for the probability generating function G(u,t) of the process N(t), t^O, the 
following representation holds for u e [0, 1] 



(2.40) 



e" st G(u,t)dt = w * 
Jo fc=0 Jo 



Pr{iV(t) = fc}di 



e~ flt Pr{ min X k > 1 - «}di 

Jo Ojjfcs;JV(t) 



fc t j n-=o(A+^) 

where Xk, k ^ I, are i.id. random variables uniform in [0, 1]. 
TTie representation of the probability generating function as 



(2.41) 



G(u, t) = Pr{ min X k > 1 - it}, 

0«fcsSiV(t) 



follows the same lines of the time and space fractional Poisson processes described in Orsingher and Polito 
[9]. In (2.40), the driving process is the state- dependent Poisson process. 

3. Alternative forms of the state-dependent Poisson process 

We construct now a point process with independent but not i.d. inter-arrival times. In particular the 
waiting time Uk between the fcth and (k + l)th arrival is distributed with p.d.f. 

(3.1) f Uk (t) = Xf^E^-Xt"*), t > 

Let us now call N(t) , t ^ 0, such process and we have the following theorem. 
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Theorem 3.1. The state probabilities of the process N(t), t ^ 0, are governed by the integral equation 

(3.2) Pk (t) -p fc (0) = -\r k Pk {t) + xr^pk-^t), t > 0, v k e (0, 1], 

where I v is the fractional integral in the sense of Riemann-Liouville (see (1.9)). Moreover their Laplace 
transforms are given by 



(3.3) 



rOO 

e- st Pr{Af(t) 
Jo 



= k}dt = A 



n u^ 3 + A ) 



Proof. First, we observe that the Laplace transform of the state probabilities, can be directly calculated 
by using the definition of the process A/"(t) 



(3.4) 



r 

Jo 



; Pr{A/"(f) = k}dt 



r 00 r rt 

e - st dt Pr([/ 
Jo L Jo 

i r 

s Jo 



U k -i e 



dy) - f Pr(J7 
Jo 



+ ■ ■ ■ + Uk e dy) 



e-" v [Pr((7 + ■ ■ ■ + £4-i e dy) - Pr(U + ■ ■ ■ + U k e dy)] 



ik+1 



1 A fc (A + s^)-A fc+1 



A ft 



We notice that, unfortunately, it does not coincide with (2.8). 

Hence we have two distinct processes that can be matched only by taking the assumption that v k = v 
for each k = 0, 1 . . . (in other words in the time-fractional Poisson case). 

We can also find in explicit way the integral equation governing the probabilities p k (t) = Pr{A/"(t) = k}. 
We start from the ordinary difference-differential equation, governing the Poisson process 

dpk , 



(3.5) 

with initial conditions 



dt 



(t) = -Apjt(t) + Ap fc _i(t), 



1 k = 0, 
k > 1. 



(3.6) Pk (0) 
By integration with respect to t, we have the equivalent integral equation 
(3.7) 



Ph(t) ~Pk(0) = -A Pk(s)ds + A p fe _i(s)ds 
Jo Jo 



In order to obtain a fractional generalization of the last equation, we replace the first-order integral in 
the right hand side of (3.7), with state-dependent fractional integrals, i.e. 

(3.8) Pk (t) - Pk (0) = -\r k Pk (t) + Xrx-^pk-^t), t>0, J/ fe e(0,l], k^O, 
where V is the fractional integral in the sense of Riemann-Liouville. For k = 0, we have 

(3.9) p (t) - 1 = -A/> (t), 

whose solution is simply given by po(t) = -E„ 0l i(— At" ). With k = 1, we obtain 

(3.10) Pl (t) = -\I vl Pl (t) + AI> (t), 
whose Laplace transform, after simple calculations, is given by 

As'" 1 " 1 



(3.H) 



Pi(t) 



(s"o + A)(s' y i + A) ' 

and coincides with (3.4) in the case k = 1. Then, it is immediate to prove that, for any order k J; 1, the 
Laplace transform of p k (t), is given by (3.4). This proves that (3.8) is the governing equation for AT(t), 
t ^ 0, as claimed. □ 
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In order to highlight the relation between the two processes N(t), t ^ 0, and J\f(t), t ^ 0, we can write 
the following 

n,=o(* 3 + A ) 

Therefore if (y^ — Vq) > for a fixed k we have that 

(3.13) p?(t) = — , f (t-j)^'- 1 ^) = fc}dj/ 

1 [y h — v Q ) Jo 

= r k ~ v ° Pr{Af(t) = fc}, 0, 

where I" k ~"° is the Riemann-Liouville fractional integral. Note that since the Riemann-Liouville fractional 
derivative (that we indicate here with D a ) is the left-inverse operator to the Riemann-Liouville fractional 
integral we also obtain the related relation 

(3.14) D" k -"«p" k k (t) = Pr{Af(t) = k}, t>0,(v k - v Q ) > 0. 
Conversely we can write 



(3.15) e~ st Pr{JV"(£) = k}dt 

Jo 



and thus if {vq — Vk) > 0, for a fixed k, we obtain that 

(3.16) Pr{AA(t) = k} = — -i s f (i - y^-^-Vk" (v) dy 

1 (^o — Vk ) J 

= r°- l,k Pl k {t), t>o, 

and that 

(3.17) D v "- Vk Pi{jV(t) = k}= p" k k (t), t ^ 0, (u - v k ) > 0. 
Finally, we have the following relation between the state probabilities of the two processes 

'f'-"»Pr(Af(!) = f:), v k >v , 
D v °~ Vh Pr{AA(t) = k}, v k < u . 



(3.18) p?{t) 



In order to deepen the meaning of this relation, we consider as an example the relation between p" 1 (t) and 
Pr{JV"(i) = 1}. 

By inverting the Laplace transform (3.4), we obtain that 

QO m / \ 

(3.19) Pr{^(t) = 1} = 2 (-ir\ m+1 £ Hr^'"-^ + ^'- +I 'o^ iii;o(m _ r)+1/ir+uo+1 (-A 2 ^ + ^) > 

m=0 r=0 \ ) 

by calculation similar to those given above for p^it). Recalling that (Mathai and Haubold [7], page 123) 

(3.20) rif^E^iat )] = x^^Efia+^atP), 
and assuming, for example v\ > uo, we find that 

(3.21) r 1 "" Pr{jV(i) = 1} 

00 m / \ 

— V I T\m\m+1 V I m I jvi-va | f »o(m-r)+«ir+i/o nra+1 / \ 2 ,u +"l \\ 
~ Zj {~ l > A Zj r I 1 I* ^ + I , 1 , l ,o(m-r) + t - 1 r + I , + l(,- A * ) j 

m = r = \ / 

00 m / \ 

— V (* ism\m+l V" 1 I m |^o(m-r-)+«i»-+''i nm+1 / >2,i/ +i/i _ ^ /,-. 

— ZjK-t) A Zj r * ^ VO + vi,vo(™-r) + vir + V1 +l(- A * J ~ Pi W> 
m = r = \ / 

as expected. 

Moreover, we observe that, since Po°(i) = Pr{7V(i) = 0} = -E„ ,i(— At" ), we have 

00 O0 

(3.22) £ pl k (t) = £ Pr{Af(t) = k} = 1 — E VoA {-\?°). 

k=l k=l 

In view of (3.18), this implies that 

(3.23) Y I Pr{Af(t) = k}=Y I pl k (t)= ^ I v "- U ° Pr{JV(t) = k] + £ D" ^ Pr{A/"(t) = fc}. 

fc — 1 fc — 1 fc: (-'fc>^o l/ k<VQ 
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The second process we construct here, denoted by N(t), t ^ 0, is given by the following generalization 
of the Poisson process, whose probabilities are given by 

(At) j 1 

r(i/,i+i) e v . i(At) 

( 3 - 24 ) Pr{7V(t) = j} = ^ +00 (xt)i " I ■ i>0, 

2jj=o r(^j+i) B„ f ,i(A0 

where A > 0, < Vj < 1. We can treat it as a generalized Poisson process with state dependent 
probabilities. Indeed we notice that, if Vj = 1, for all j, we have 

(3.25) Pr{%) = j) = - Jgg - = ^e"" = Pr{N(t) = j}, 

that is the state probability of the homogeneous Poisson process. 

A similar construction was adopted in Beghin and Orsingher [3]. We notice that an analogous general- 
ization was used by Sixdeniers et al. [13] in quantum mechanics, in relation to Mittag-Leffler type coherent 
states. We now recall from Balakrishnan and Kozubowski [1] that the distribution (3.24) can be regarded 
as a weighted Poisson sum. Indeed we notice that 



(3.26) 

Hence we have 



(Xty 1 _ f(^+i)Pr{N(i)=j} 



rfa j + 1) E v , A {\t) T^fhlJ Pr{N(t) = k} ■ 



Pr{N(t)=j} 



(3.27) Pr{N(t)=j} = + ] 



w 3 ' ■ Pr{N(t)=j} 
Zjj =0 " 



JJ=U r ( ^ + i) P r {N(t) = fc} 

The probability generating function of (3.24) is given by 

GO 

(3.28) G(u,t) = u k Pr{N{t) = k} 

k = 

T+GO (\ut) k 



s 



fc=0 r^fc+l) £„ fci i(At) 
(At)* l 

r(v k k+x) B„. ,i(At) 



V+°° (At) fc 1 

Zjfc=o 

In the case Vj = v, for all j ^ 0, we have 

(3.29) G(u,t) = ^ ^ }» = 

Zjfc=0 r(i/fc + l) E„ a (At) 

that coincides with the equation (4.4) of Beghin and Orsingher [3]. 

By means of the generating function we can also find the explicit form of the mean value of the distribution 
(3.24), i.e. 

\j-V + ° k(\t) k - 1 1 \tY + CC ' (At) fc 1 

(3 30) EN(t) = k=0 r ('^ fc+1 ) g ^,i( At > = ^ k=0 ^+i r (^+i fc +^+i) ^ t+1 ,i(At) 

''''' ^ ' (Ai) fc l V+ 00 (At) fc 



Zjft=0 rfi/i.fc+1) E„, ifAtl ^ 



rc^fc+i) B^,i(At) ^fe=o rfofc+i) s„ fe ,i(At) 

such that, in the limiting case Vk = v for all k, we recover the case considered in Beghin and Orsingher [3] 
and in Beghin and Macci [2], i.e. 

We now consider a sequence of a random number of non-negative i.i.d. random variables with distribution 
F(fS) = Pr(X, < /3), i ^ 1 and represented by N(t). The distribution of the maximum and minimum of 
this sequence is given by 

Y+co (AF(/3)t) fc 

(3.32) Pr{ma X (x 1 ,-,X %) )<fl = 



vn + oo 
2jfc = 



r> fc fc + l) g„ fc ,i(At) 

(Ai) fc 1 
r> fe k+l) i3„,,i(At) 
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+<* (\[1-F(f3)]t) k 



(3.33) Pr{min (x 1; ■ ■ ■ ,% (f) ) > /3} = — 



r(v k k+l) £„, i(At) 



yi+oo (Xt) k 1 
2-ik=0 rO fc fc+l) E„ ,i(Ai) 



In the case = = 1, for all k, we recover the distribution of the maximum and minimum of the 
homogeneous Poisson process. 

4. State dependent fractional pure birth processes 

In this section we consider a different point process which can be generalized in a state-dependent sense 
as we have done for the fractional Poisson process. We thus analyze a state-dependent fractional pure birth 
process (see Orsingher and Polito [8] for the fractional case with constant order), where the probabilities 
are governed by the following equations 



f ^Pk k (*) = -AfcPfc* W + Xk-ipl^it), k > 1, t > 0, v k 6 (0, 1], 



(4.1) 



dt* 

Pl k (0) = 



1 k = 1, 

k =s 2, 



As in the Section 2 the Laplace transform of the solution to (4.1) can be found rather easily. This is 
done in the following proposition. 

Proposition 4.1. The Laplace transform of the solution to the state- dependent fractional pure-birth process 
(4.1) reads 



'(«) = \ + J e- st pl»{t)At= \Y\\ 



(4.2) p k k '^- 
where the fractional derivative appearing m (4.1) is in the sense of Dzherbashyan-Caputo. 



Proof. We can solve equation (4.1) by means of an iterative procedure, as follows. The equation related 
to k = 1 

(43) f d wK 1 W = -Ai P r i W, t>0,i*e(0,l], 

1^(0) = 1, 

has solution p" 1 (t) = E vi ,i(—\t vi ). For k = 2, the equation 



(4.4) 



= -\2P?(t) + Al P ?(t), t>0,v 2 e (0, 1], 

[p?(0) = 0, 
has solution with Laplace transform 



(4.5) P?(s)=\ e- st p?(t)At= ^ — 1 

Jo Ai + s"i A 2 + 

whose inverse is given by (see (2.9)) 



I 

oo m / \ 

(4.6) p-w = ^ (-i)'" ^ ™ A; +1 Ar- r ^-' + - + x + 4„( ro - r)+ , rrt+1 (-AiA 2 r i+ " 2 ). 

m = r=0 V, / 

By iterating this procedure, we arrive immediately at 

(4.7) p v * (s) = J"*" e-'Yj (t)dt = iff A 



' ' n-=i(^+Ai)' 

as claimed. □ 

By recalling (2.28), we obtain the explicit expression of the state probabilities p^. h (t), k > 1, t > 0, as 

(4.8) pj*(t) = £„!,!(- Ait" 1 ) * | W^Ev^A-Xjt"'). 

We now consider the state dependent linear birth process, denoted by Nu n (t), t ^ 0. This means that 
we take A& = Afc in (4.1). We have the following 
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Theorem 4.2. Let us consider the state dependent linear birth process Nu n (t), t ^ 0, governed by 

^Pl k (t) = -AfcpJ»(t) + \{k - l)p£?{t), k>l,t>0,v k e (0, 1], 
( 4 -9) \^ k <^ fl *=L 



Pfe fc (0) = 
the following relation holds 
(4.10) 



k > 2, 



fe=l J=l \ J / 



m—j + 1 



Proof. In order to find explicit relations for the moments of the distribution Nu n (t), we multiply both 
sides of equation (4.9) by fc m and sum over all the states, obtaining 

00 CO CO 

( 4 - n ) £ * m ;L-#(*) = ~ A S ^W*) + a s - i) ft ^(t) 



-A ^ fc m+ V (t) + A ^ fc(fc + 1)>^ (t) 



fe=l fc=l 

oo oo m 



-A ^ JT +1 P ? W + A^ . fc— (t) 

fe=l fc=l j=0 \3 ) 

m / \ 00 m / 

i=l V- 7 / fc=l J=l V 



□ 



Remark 4.3. We can consider in a explicit way the relations involving first and second moments. For 
example, if we multiply (4.9) for k and sum over all the states, we obtain that 

CO CO CO 

(4-12) £ k±-p? (t) = -A £ fc 2 P ^ (t) + A 2 k(k - l)^! 1 (t) 

fc=l fc=l fc=l 

OO 

= A^fep^(f)=AEiVi <n (t). 

fc=0 

In the same way, for the second moment, we multiply (4.9) for k 2 , obtaining 

CO n„ fc CO CO 

(4-13) ^ k 2 ^P?(t) = -A X fc 3 p?(*) + A £ fc 2 (fc - lJp^iH*) 

fc=i fe=i fe=i 

= Af]fcp^( i )+2Af;fe 2 p^( i ) 

fc=0 fe=0 

= XEN Hn (t) + 2XE(N Un f(t). 
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